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We present an experimental and theoretical study exploring surface effects on the evolution of 
the metal-insulator transition in the model Mott-Hubbard compound Cr-doped V 2 O 3 . We find a 
microscopic domain formation that is clearly affected by the surface crystallographic orientation. 

Using scanning photoelectron microscopy and X-ray diffraction, we find that surface defects act as 
nucleation centers for the formation of domains at the temperature-induced isostructural transition 
and favor the formation of microscopic metallic regions. A density functional theory plus dynamical 
mean field theory study of different surface terminations shows that the surface reconstruction 
with excess vanadyl cations leads to doped, and hence more metallic surface states, explaining our 
experimental observations. 

PACS numbers: 73.20.At; 


Metal to insulator transitions (MIT) are among the 
most remarkable macroscopic effects of electronic corre¬ 
lations in condensed matter. After many experimental 
and theoretical studies, it has been possible to under¬ 
stand the crucial role played by the lattice, which can 
stabilize electronic instabilities and guide the evolution 
of correlation-driven phenomena such as Mott-Hubbard 
transitions [T]. However these studies have concentrated 
on bulk properties: the surface behavior is rarely dis¬ 
cussed and is indeed a more complicated problem. At 
the surface, the atomic coordination number and screen¬ 
ing change and this has an effect on electronic correla¬ 
tions, but other factors such as surface reconstruction 
and lattice defects can also affect the MIT. 

Cr-doped V2O3 is the prototype Mott-Hubbard ma¬ 
terial US], presenting a correlation-induced MIT with¬ 
out symmetry breaking. Strong electronic correlation 
split the noninteracting bands into interacting upper and 
lower Hubbard bands. In the metallic phase a strongly 
renormalized quasiparticle peak remains at the Fermi 
level which is reminiscent of the uncorrelated band struc¬ 
ture. The phase diagram of V2O3 consists of three 
phases: paramagnetic insulator (PI), paramagnetic metal 
(PM), and antiferromagnetic insulator (AFI) [ 2 ] (see Fig. 
1 a)). The Mott transition takes place between the PI 
and PM phases, and it can be induced by increasing 
pressure starting from the PI phase, as well as by de¬ 
creasing temperature for doping levels around 1 . 1 % Cr 
concentration. The orbital degrees of freedom in V2O3 
have to be taken into account in order to understand 
the MIT: the low-lying orbitals are the t2^ orbitals made 
up of singly degenerate ai^ and doubly degenerate e^ 
orbitals in the crystal field of the corundum lattice struc¬ 
ture. At the MIT the ai^ bands are shifted up in energy 


and the e^ bands split into two Hubbard bands, with 
the Hunds exchange leading to a local spin alignment 
throughout the transition ISHH]. Although the surface 
of V2O3 has been well studied for the pure compound 
mm, and extensive experimental [12] and theoretical 
investigations m have been carried out on the surface 
termination of V2O3, very little is known about the ef¬ 
fects of the surface on the Mott transition. In this let¬ 
ter, we explore specific surface effects on the evolution of 
the Mott-transition in (Vi_a,Gra,)203 both experimen¬ 
tally (using X-ray diffraction (XRD) and scanning pho¬ 
toelectron microscopy (SPEM)) and theoretically by den¬ 
sity functional theory plus dynamical mean field theory 
(DFT+DMFT) gOHil] 

The Mott transition is a first order transition there¬ 
fore domain coexistence can occur and a large hysteresis 
is possible. Fig lb) shows XRD measurements on a 
monocrystal done during the temperature driven phase 
transition m- A large hysteresis is seen comparable to 
earlier results [ 3 ] . However the transition happens in less 
than 2 K at 244 K whereas in powder XRD the transition 
lasts 60 K. The major difference between a monocrystal 
and a powder is the surface to bulk ratio. The MIT at 
the surface might be greatly effected by defects or surface 
reconstruction. In order to understand this discrepancy 
we probe the surface using SPEM. 

Method - Experimental. SPEM experiments were per¬ 
formed as a function of temperature on the Escami- 
croscopy beamline at Elettra synchrotron, using photons 
of 400 eV [To]. SPEM uses a direct approach to photo¬ 
electron spectromicroscopy which is the use of a focused 
photon beam, down to 150 nm in diameter, to illumi¬ 
nate the sample. Photoelectrons are collected with a 
hemispherical electron analyzer and detected by a 48- 
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FIG. 1: a) Phase diagram for (Vi-ccCrcc )2 03 . b) Ratio 
PM/(PM+PI) versus temperature probed by X-ray diffrac¬ 
tion on a monocrystal. c) Photoemission spectra for 
(Vi_xCraj )203 at x=0.011, 300 K and 200 K. The increase 
of spectral weight near the Fermi level indicates the phase 
transition from insulator to metal when decreasing the tem¬ 
perature. and Ilhb represent respectively the intensity 
near Ef and of the LHB. The ratio between the two provides 
a measure for the metallicity of the probed surface. 


channels electron detector. SPEM can operate in two 
modes: i) XPS spectroscopy from a sub-micron spot; ii) 
imaging where the sample surface is mapped within a se¬ 
lected kinetic energy. In imaging mode, all 48 channels 
are recorded for every pixel of the image. In our exper¬ 
iments a 2 eV range near the Fermi level was chosen to 
observe the MIT, as seen in Fig. 

As it has already been shown in previous SPEM 
measurements m, coexisting PM and PI phases in 
(Vi_a^Cra,)2 03 can be unambiguously distinguished by 
the photoemission signal from the outer valence elec¬ 
tronic states, but the physical origin remained myste¬ 
rious. With the present study we are able to solve this 
puzzle. The improved experimental conditions allow us 
to show a clear correlation between the shape and po¬ 
sition of domains at the metal-insulator transition, and 
the fact that surface structural defects actually favor the 
formation of metallic domains. This was possible thanks 
to: i) the use of higher energy photons, that make it 
possible to have a clearer spectroscopic contrast [H]; ii) 
an improved spatial resolution (150 nm); and iii) a di¬ 
rect comparative study between different surfaces of the 
same material. In particular, high quality single crys¬ 
tals from Purdue University were either cleaved in the 
(001) plane or the (102) plane, keeping always the [001] 
direction pointing towards the electron analyzer which 
was shown to give the strongest quasiparticle peak in 
the photoemission signal [ 19 ] (see Fig. 2 top). The sam¬ 
ples were cleaved and measured under UHV conditions 
(2x10“^^ mbar or better) to avoid surface contamination 




FIG. 2: Top: Sample orientation versus incoming photons 
and analyzer Left: Sample cleaved along the (102) direction 
with analyzer in (001) direction, a): Surface morphology 
image taken at a photoelectron energy corresponding to the 
V3p core level, and at temperature T = 245 K. b) contrast 
IefIIlhb between the PM and PI phases for 245 K and 
c) for 243 K. No coexisting domains could be detected dur¬ 
ing the phase transition. Right: Sample fractured along the 
(001) direction with analyzer in (001) direction, d) Surface 
morphology image. The following images show the contrast 
between the PM and PI phases: e) T=267 K, f) T=257 K 
and g) T= 242 K. 


(the carbon photoemission peak was regularly checked in 
the course of the experiment). In order to obtain genuine 
two-dimensional maps of the PM vs PI concentration at 
the surface, the photoemission images were corrected for 
topographic effects [20]. As a measure for the metallic¬ 
ity of the sample, we take the ratio between the intensity 
near Ep 0 -Ef) near the Lower Hubbard Band (LHB) 
0-lhb), see Fig. 0 The images obtained in this way are 
real-space ’’metallicity maps” of the sample surface m- 
SPEM results Fig. shows a region on the ( 102 ) 
cleaved surface at different temperatures while cooling 
down the sample. Since the (102) plane is the natural 
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cleaving direction, it can produce large flat surfaces with 
few defects. In particular, the probed region does not 
present any significant defects detectable with our spatial 
resolution (150 nm): this is evident from Fig. |^a), where 
the image has not been corrected for surface morphology 
and which cannot reveal the presence of any topographic 
defect at the surface. By observing the metallicity im¬ 
ages we found that no MIT could be seen until 245 K 
(Fig. i b)), while the whole surface turned metallic at 
243 K (Fig. I^c)). This indicates that the transition was 
too sudden for our experimental conditions to detect the 
presence of domains, the limiting factor being the rela¬ 
tively long acquisition time for each photoemission image 
and/or the size of the domains. Hence, we can conclude 
that, within the limits of our spatial resolution, the entire 
probed region turns metallic in less than 2 K, similar to 
the bulk probed by our XRD experiment. Therefore the 
( 102 ) plane behaves like the bulk material. 

On the contrary, when cleaving along the (001) direc¬ 
tion, we find coexisting domains on the surfaces over a 
wide temperature range, starting as soon as 267 K. We 
start by analyzing an area where no structural defects 
could be detected with our spatial resolution, as revealed 
by the uncorrected image in Fig. [^d): the first metallic 
domains appear starting at 267 K (Fig. [^e)), and they 
present a triangular shape, revealing a clear correlation 
with the hexagonal symmetry of the ( 001 ) plane (the a 
and b crystallographic directions are shown to allow a di¬ 
rect comparison): this strongly suggests that the borders 
of the metallic domain correspond to surface steps. The 
rest of the area is still insulating and starts its transfor¬ 
mation around 257 K (Fig. [^f)). Eventually at 242 K 
(Fig. [^g)), the transition is almost complete, similarly 
to the ( 102 ) plane. On the (001) surface, the sample also 
has regions with a significant amount of large defects as 
represented in Fig. Here large metallic domains ap¬ 
pear between the cracks. Let us now focus our attention 
on smaller defects, most likely corresponding to cleavage 
steps, such as the one found in the ”S” area shown in de¬ 
tail in Fig. [^c): the metallic domains appear to follow 
the cleavage step direction. Overall, metallic domains are 
present for the (001) surface well above 243 K, which is 
the transition temperature of the ( 102 ) plane. 

Our experimental findings confirm that structural de¬ 
fects can act as nucleating centers for the insulator-to- 
metal transition, thus guiding the evolution of the for¬ 
mation of domains during the coexistence [16]. Surpris¬ 
ingly, they also clearly indicate that structural defects 
at the surface of (Vi_a;Cra ,)203 favor the formation of a 
metallic phase over an insulating one. 

Theory. One might envisage the cleavage steps as an 
extra surface so that surface effects should be enhanced. 
However the usual suspects of electronic correlations at 
(enhanced) surfaces, i.e., an enhanced U because of the 
reduced screening and the removed hopping perpendicu¬ 
lar to the surface, suggest the surface and cleavage steps 




FIG. 3: Sample cut along the (001) plane at 267 K. a) Sur¬ 
face morphology image with a crack marked by the dashed 
arrows and a smaller defect marked by the S box. b) The 
corresponding phase contrast images show the apparition of 
metallic domains starting at 267 K. The large domains are 
separated by cracks in the sample, c) Close up image of a 
cleavage step, which presents a more metallic behavior, d) 
Photoemission spectra from the Pi and P 2 points. 


to be more insulating. We find the contrary; they are 
more metallic. 

With the simple explanations at failure, we need to in¬ 
spect the surface more thoroughly, including the proper 
surface reconstruction. DFT calculations [39] indicate, 
depending on the oxygen partial pressure, a VO surface 
termination with various excess of vanadyl or a O 3 ter¬ 
mination. The former has also been identified in some 
experiments HIHIllIl]. Fig. 4 shows the actual super¬ 
cells with VO and O 3 termination that we consider in 
our calcualtion. One can envisage the surface as consist¬ 
ing of several layers of the stoichiometric (V-O 3 -V) unit 
cell, i.e., (V-O 3 -V) (V-O 3 -V) ••• (V-O 3 -V) plus an ex¬ 
tra V or (O 3 -V) layer for the VO and O 3 termination, 
respectively, see Fig. 4 (upper panel). 

Method - Theory. We perform a full relaxation of the 
VO and O 3 terminated V 2 O 3 surfaces in the corundum 
structure including a 12 A thick vacuum layer, using 
the Vienna ab-initio simulation package (VASP) [30] [3T]. 
with GGA-PBE functional [32]. In the case of the O 3 ter¬ 
mination, the large polarity of the surface is compensated 
by a surface reconstruction where one subsurface V atom 
moves from the second layer to the first (surface) layer 
so that (O 3 -V) (V-O 3 -V) (V- 03 -V)--- actually becomes 
(O 3 -V 3 -O 3 ) (V- 03 -V)''-, the surface reconstruction for 
the VO termination is less dramatic. 

After atomic relaxation we Wien2Wannier-project [34] 
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the corresponding Wien2K [33] bandstructure and sup¬ 
plement it by a local Kanamori interaction with intra¬ 
orbital U = 5.5 eV, inter-orbital U' =4.1 eV and Hand’s 
exchange J = 0.7 eV. The interaction parameters are 
taken a bit larger than in the literature iSHaEI] first be¬ 
cause of the reduced screening at the surface and second 
to account for the more insulating nature at the exper¬ 
imental 1.1% Cr-doping, which is too small to take into 
account in the supercell of our calculations. The resulting 
Hamiltomian is then solved by DMFT at room tempera¬ 
ture, 300 K, using continuous-time quantum Monte Carlo 
(CTQMC) simulations in the hybridization expansions 
[36] (w2dynamic code [37|) and the maximun entropy 
method [38] for an analytic continuation of the spectra. 

DFT+DMFT results. Fig. 4 (lower panels) shows the 
DFT-fDMFT spectra for VO (left) and O 3 termination 
(right), resolved for the aig and orbitals and the dif¬ 
ferent layers. The VO terminated surface is insulating at 
this interaction strength and temperature. Compared to 
a 3d^ electronic configuration for all V atoms, the vanadyl 
termination adds 1 extra O to each (V-O 3 -V) unit cell in 
the layer, i.e., 1 hole per V in the layers. According to 
our DFT+DMFT results, this hole is however bound to 
the surface layer where the single V in the surface layer 
is in a V^+ or 3d^ configuration and to the second layer 
where one of the two V atoms is V^+. The other V atom 
in the second layer and all other V in the further sub¬ 
surface layers are in a V^+ or 3d^ configuration. This 
charge disproportionation explains why, despite of the 
doping, the surface may remain insulating. (Note, Fig. 
4 show the layer-averaged spectrum, the supplemental 
material resolves it for the two inequivalent sites of the 
second layer.) If we reduce the interaction strength, the 
VO terminated surface becomes metallic. Indeed the VO 
terminated surface is more metallic than bulk V 2 O 3 , i.e., 
it stays metallic up to larger values of the Coulomb in¬ 
teraction (see Supplemental material) 

The O 3 termination (Fig. 4 right) is, on the other hand, 
already metallic at the same Coulomb interaction and 
temperature. The most metallic layer is the surface layer; 
and the width of the central quasiparticle peak shrinks 
from layer-to-layer. The fourth, central layer is already 
close to the bulk result (see Supplemental Material). The 
reason why the surface layer is more metallic, despite the 
reduced hopping, is the even larger hole-doping due to 
the O 3 termination. The (O 3 -V 2 ) slab adds to the stoi¬ 
chiometric (V-O 3 -V) layers in the case of the O 3 termi¬ 
nation 3 holes per (V-O 3 -V) unit cell in the layers (or 
1.5 holes per V atom). These 1.5 holes are now however 
distributed to three surface layers in the DFT+DMFT 
charge distribution: 1.2 electrons per V (first layer), 1.5 
(second layer), 1.9 (third layer), 2.00 (i.e., the bulk value 
for the fourth layer). That is for the O 3 termination 
Vanadium is neither nor V^+ but in between. Con¬ 
sequently, the system is more itinerant; the surface layers 
of V 2 O 3 are metallic. 
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FIG. 4: Top: VO termination (left) and Supercell O3 termi¬ 
nation (right). Below the figures, we indicate how the super¬ 
cell is made up from the stoichiometric (V-O3-V) supercell 
plus the surface termination. For the O3 termination, one 
V moves from the second to the first layer because of the 
surface reconstruction. Bottom: Layer- and orbital-resolved 
spectral function A{uj) (a)-(b) for the VO termination, (c)- 
(d) for the (O3) termination. The more excess oxygen at the 
surface, the more metallic are the surface layers. 


Hence we can conclude that excess oxygen makes the 
(reconstructed) surface more metallic due to hole-doping. 
The additional surface at the steps indicates an extra 
surface doping so that V 2 O 3 will be even more metallic 
at such edges. This is akin to our O 3 termination vs. 
VO termination. As we have seen, this larger amount of 
doping has not a very big effect, but it shifts the criti¬ 
cal Coulomb interaction for the Mott-Hubbard transition 
somewhat to the left of the phase diagram (inset Fig. [^. 
That means the transition temperature increases as we 
found experimentally. It also explains our experiment. 
Fig. i the cleave step is more metallic because of an 
accumulation of excess oxygen at the corner. 

Conclusions. We observed experimentally that at the 
PI/PM Mott transition in (Vi_a^Cra ,)203 metallic do¬ 
mains appear at higher temperatures than the bulk tran¬ 
sition. Their evolution is determined by the surface crys¬ 
tallographic direction and along the cleaving steps. Our 
DFT+DMFT theoretical calculations show that a sur- 
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face reconstruction with an excess of oxygen favors the 
formation of a metallic phase. Against common wisdom, 
sufaces can hence be more metallic than the bulk, and 
surface steps even more so. This effect observed here 
for (Vi_a;Cra ,)203 can be of general interest for surfaces 
of strongly correlated oxides, oxide heterostructures, and 
nanostructure. It is also relevant when comparing surface 
and bulk sensitive experiments. 
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